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ABSTRACT 

A three dimensional Smoothed Particle Hydrodynamical (SPH) model has been extended 
to study the radiative driven implosion effect of massive stars on the dynamical evolutions 
of surrounding molecular clouds. The new elements in the upgraded code are the inclusion 
of Lyman continuum in the incident radiation flux and the treatment of hydrogen ionisation 
process; introducing ionisation heating & recombination cooling effects; and adding a proper 
description of the magnetic and turbulent pressures to the internal pressure of the molecular 
cloud. This extended code provides a realistic model to trace not only the dynamical evolution 
of a molecular cloud, but also can be used to model the kinematics of the ionisation & shock 
fronts and the photo-evaporating gas surrounding the molecular cloud, which the previous 
code is unable to deal with. 

The application of this newly developed model to the structure of the middle Eagle Neb- 
ula finger suggests that the shock induced by the ionising radiation at the front side of the head 
precedes an ionisation front moving towards the center of the core, and that the core at the fin- 
gertip is at transition stage evolving toward a state of induced star formation. The dynamical 
evolution of the velocity field of the simulated cloud structure is discussed to illustrate the role 
of the self-gravity and the different cloud morphologies which appear at different stages in the 
evolutionary process of the cloud. The motion of the ionisation front and the evaporating gas 
are also investigated. The modelled gas evaporation rate is consistent with that of current other 
models and the density, temperature and chemical profiles are agreement with the observed 
values. 

The relative lifetimes of different simulated cloud morphologies suggests a possible an- 
swer to the question of why more bright-rimmed clouds are observed to possess a flat-core 
than an elongated-core morphology. 

Key words: star: formation - ISM: evolution - ISM: HII regions - ISM: kinematics and 
dynamics - radiative transfer 



1 INTRODUCTION 

The Eagle Nebula (Ml 6) presents an excellent laboratory in which 
to study the influence on star formation due to the presence of 
nearby massive stars, which can lead to radiative implosion of sur- 
rounding molecular cores and the initiation of secondary star for- 
mation (Bertoldi 1989; Lefloch & Lazareff 1994). The intriguing 
structures in the heads of the Elephant Trunks in the Eagle Nebula 
have inspired many studies. A comprehensive survey conducted by 
iHester et all jl99a) using the HST Wide Field and Planetary Cam- 
era (WFPC2) resolved a number of Evaporating Gaseous Globules 
(EGGs) at the tip of finger-like features protr uding from column s 
of cold gas and dust of the Eagle Nebule iHester et alJll99^ . 
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More investigations on the physical characteristics and the future 
evolu tion of the head s t ructures have also been carried out since 
then jWhite et alJll999t IWilliams et al]l2()0lt Ipukuda et alJl2002t 
ISugitani et all2002l:lThompson et alj2002l:IUrauhart et ai2003ir ^ 

In one of our early papers jwWteetn][l99i), we reported 
molecular line, millimetre/submillimetre continuum, and mid-IR 
observations of this region. The steady state of the head structure of 
the middle finger of the Eagle Nebula was estimated using a SPH 
model which included a simple treatment of an isotropic FUV radi- 
ation field (6.5 < hv < 13.6 eV). The comparision of the estimated 
steady state of the head structure of the middle finger with the ob- 
served data suggested that the FUV ionisation induced shock might 
have not passed the core at the middle fingertip, so that the head 
structure of the fingertip could not be taken as the remains of ra- 
diatively driven implosion. The observations revealed that the core 
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structure shows the characteristics similar to those expected in the 
earliest stages of protostellar formation, therefore it was suspected 
that the head structure of the middle finger in the Eagle Nebula 
might evolve from a large and dense structure which existed prior 
to the expansion of the HII region. 

Recently there have been other observations that have sup- 
ported the idea of triggered protostellar formation at the tip 
of the middle finger, following the interaction of the natal 
molecular cloud with the ionising radiation from the nearby 
OB stars jpukuda et all I200I iMcsaughrean & AnderserJ booa 
ISugitani'etanT200^ r iTliomDsoi^^TT2W The question is 
whether or not UV radiation from the nearby massive stars could 
trigger the collapse of the head structures in Eagle Nebula and lead 
to the next generation of star formation. The challenge is to develop 
theoretical models that provide a reliable description of the dynam- 
ical evolution of these clumpy structures and which can trace back 
the history of the clumps, accommodate the current observations, 
and predict the future evolution of these fingertip structures. 

A precise description of the genesis and the future evolution 
of the head structures such as those of the Eagle Nebula fingers re- 
quires a more rigorous dynamical modelling which should include 
Lyman continuum in the external radiation field, since hydrogen 
ionisation ought to play a dominant role in the dynamical evolution 
of the molecular cloud structures like those in Eagle Nebula. 

A substantial amount of theoretical work has already been 
carried out to investigate the evolution of dense, gaseous clumps 
which are in the vicinity of massive OB stars. T hese include the 
Radiatively Driven Implosion (RDI) mechanism jBertoldilfTgsl 
which was developed using an approximate analytical solution for 
the evolution of a spherical symmetric neutral cloud subject to the 
radiation o f a newly formed nearby star; and the cometary glob- 
ule model <Lefloch & Lazarefjll994 \\995.\ which describes the 
dynamical evolution of a neutral globule illuminated by the ion- 
ising radiation of OB stars based on 2-D hydrodynamical simula- 
tions. Although these 2-D models successfully reproduced some of 
the observed characteristics of objects influenced by nea rby mas- 
sive s tars, they did not adequately include self-gravity iBertoldi| 
Il989^ , nor self-consistently treat the chemistry and thermal evo- 
lution iLefloch & La zarefl:^1997.'l995'l. Consequently they are un- 
able to deal with the details of triggered star formation by UV radi- 
atioii; 

IWilliams et alj i200lh presented a 2-D (i.e., cylindrical sym- 
metry of the cloud structure is used) hydrodynamical simulation of 
the evolution of photoionised clumps with similar characteristics to 
those found in the Eagle Nebula. The simulated results over a vari- 
ety of initia l conditions propos ed a different evolutionary scenario 
from that of|wWteetdJ|l99i) for the head structures, suggesting 
that the head structures at the tips of the Eagle Nebula trunks may 
have already survived the propogating shock induced by the exter- 
nal UV radiation field and are in a near-equilibrium state about 10' 
years later. In this modelling, self-gravity of the molecular gas was 
also not included. 

The existing discrepancy in the descriptions of the structure 
of the heads at the finger tips and the recently reported new ob- 
servations on the finger tips' structure suggests the necessity of re- 
visiting the issue of UV radiation triggered star formation at these 
fingertips. In order to do this, a self-consistent and comprehensive 
3-D BRCs' model is required which could be used to study the UV 
radiation triggered star formation in BRCs and which more gener- 
ally can be used to describe observations in all of BRCs, proplyds 
and HII regions, to meet with the ever-increasing interest in the in- 
vestigation of the nature of BRCs and proplyds in HII regions. It 



is obviously seen from the above that this fully 3-D BRCs model 
should be able to self-consistently and simultaneously treat the self- 
gravity, thermal evolution, radiative transfer and chemistry. There- 
fore it is our intentin to develop the first self-consistent and com- 
prehensive 3-D BRCs's model based on a previous 3-D SPH code, 
which was originally develped to investigate the evolution of iso- 
lated molec ular clouds under the effect of the interstellar medium 
radiation (Nelson & Lan2er'l997). 

In the following sections in the present paper we will firstly 
present an outline for the previous 3-D SPH model and then in- 
troduce new elements neccessary for simulating the evolution of 
BRCs. The simulation results will be discussed and compared with 
observations, and the evolutionary of the simulated cloud structure 
will be investigated. A comparision of the resuluts with those based 
on the previous code is made and then we infer the future evolu- 
tion of the Eagle Nebula finger from its current state. Based on the 
good agreement between simulations and observations, a reason- 
able explanation for the observed unbalanced numbers between the 
flat-core and elongated-core BRCs is proposed. 



2 THE THREE DIMENSIONAL SPH MODELLING 

The physical properties of the head structure of the Eagle Neb- 
ula presented in our previous paper (White et al. 1999) was based 
on the Smoothed Particle Hydrodynamics (SPH) code developed 
bv Nelson & Langer 11997.) . which self-consistently treats the self- 
gravitation, chemical and thermal evolution in the molecular cloud. 
Starting with a cylindrical finger of radius 0. 1 pc with a spherical 
molecular cloud of 3 1 Mq at its head, which is under the effect of 
a FUV radiation field (with photon energy lower than 13.6 eV and 
enhanced from the +z-direction). Based on an assumed static den- 
sity distribution (i.e, under equilibrium assumption), the tempera- 
ture and chemical distributions, after the FUV radiation is swiched 
on for 10' years, showed similar characteristics to observations of 
the middle finger. However the detailed and accurate description on 
the whole dynamical evolution, especially on the genesis and future 
evolution of the observed structures was not discussed at that time, 
since the original SPH code did not include the hydrogen ionising 
influence of the Lyman continum radiation {hv > 13.6 eV), the ioni- 
sation heating, recombination cooling or the relevant radiative tran- 
fer process. These are clearly important since the effect of the UV 
radiation of nearby massive stars on their surrounding molecular 
clouds tells us that hydrogen ionisation is one of the most important 
physical processes in radiation driven implosion mo dels iBertoldi 
ll98StlLefloch & Lazareflil994IWimams et all200ll) 

Recentlv. I lKessel-Devnet & Burkerlll200(ll2003ft have devel- 
oped a three-dimensional SPH model of the radiation-driven im- 
plosion of molecular cloud which includes both self-gravitation 
and hydrogen ionisation from nearby stars. Their model provided 
the first 3D SPH code to treat ionising radiation in turbulent astro- 
physical fluid flows, opening a wide field of applications involving 
feedback processes of young massive stars on their parental clouds. 
However a self-consistent treatment for the thermal evolution was 
not included in this model and the energy evolution was approx- 
imated by a simple linear function of the ionisation ratio x, i.e., 
E = xE]Qooo + (1 - ^)Eciouii, where E is the internal energy of a 
representive particle with ionisation ratio x ( £10000 being the in- 
ternal energy of a fully ionised particle at a temperature of 10000 
K and Edoud being the internal energy of a neutral particle at a 
temperature of 10 K). The absence of a self-consistent treatment 
of the ionisation heating and recombination cooling processes in 
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their model prevents a rigorous calculation of the physical proper- 
ties of the cl oud, such as gas temperature, density, sound velocity 
and pressure j Kessel-Devn et & Burkert"2000'. 2003), which in fact 
can adequately be validated through observational data. 

We therefore extend our previous code to include a wider 
range of UV radiation energy distributions, hydrogen ionisation, 
radiative transfer, ionisation heating and recombination cooling, to 
provide a more comprehensive and self-consistent investigation on 
the effect of radiation of massive stars on their parental molecular 
clouds. 



observe the influence of the interstellar radiation on the chemical 
and thermal evolution of an isolated molecular clouds, such as Bok 
globules. For the investigation of the influence of more intensive 
UV radiation from massive stars on the evolution of its surround- 
ing objects, hydrogen ionisation heating will play a dominant role 
since the intensive Lyman continum photons available in the ra- 
diative flux, and the much higher abundance of hydrogen atoms 
than that of grain particles. Therefore it is necessary to include the 
hydrogen ionisation process and the relevant heating and cooling 
process in order to fulfill our goal. 



2.1 The main components in the previous SPH code 

A detailed description o f the previous SPH cod e can be found in 
iNelson & LangeJi 19971) and lWhite eTai] ^99^. Here we present 
a brief introduction of its main features for completeness. In the 
previous version of the code, the SPH numerical technique was 
employed to solve the following continuity, momentum and energy 
equations for a compressible fluid 







(1) 



(2) 



(3) 



(4) 
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and the chemical rate equations take the general form: 
dX, 

where ^ = ^ + v • V denotes the convective derivative, p is the 
density, v is the velocity, P is the pressure, Syisc represents the vis- 
cous forces, tl is the internal energy per unit mass, and O is the 
gravitational potential. F and A represent nonadiabatic heating and 
cooling functions respectively. 

The heating function F in thermal modelling is mainly pro- 
vided by the photoelectric emission of electrons from grains illu- 
minated by the incident FUV (6.5 < hv < 13.6 eV) component 
at the surface of the cloud. This is stronger on the surface facing 
toward the source star, and weaker on the shielded cloud surface 
since the latter is mainly from backscattered light by the surround- 
ing interstellar medium and is 0.1- 0.2 times that on the directly 
illuminated surface I Des ert et a l. 1990; Hurwitz et al. 199y, i.e., 
the FUV radiation flux incident the boundary of the cloud can be 
written as 



JFuv{r = R,e): 



2OOOG0 
0.2 * 2OOOG0, 



0° < < 90°) 
90° < < 180°) 



where Go is the standard interstellar UV field (Habina 196^) and 6 
is the azimuthal angle in spherical coordinate. Heating of the gas 
is also afi^ected by cosmic ray heating, H2 formation heating and 
gas-dust thermal exchange. 

The cooling function A is affected by CO, CI, CII and OI line 
emission. The fractional abundance of main chemical species X, in 
Equation are for CO, CI, CII, HCO+, OI, H^, HJ, OH„ CH„ 
and electrons; A", is the associated chemical reaction rate and 
n is the total number density. Detailed formula for these physical 
and ch emi cal processes can be found in papers of Nelson & La nge3 
|l993) and lWhiteetai]h999h . 

The dominant heating process included in the above is the 
photoelectric emission of electrons from dust grains by the inci- 
dent FUV (1-13.6 eV), as the code was originally developed to 



2.2 New elements 

The new elements added into the above original code are solv- 
ing for the ionising radiation transfer equation, hydrogen ionisa- 
tion heating and electron recombination cooling. We also introduce 
a realistic description of the pressure of the molecular cloud by 
including the turbulent and magnetic pressures. More detailed de- 
scriptions are contained in the following sub-sections. 



2.2.1 Ionising radiation transfer 

Although Helium ionisation was included in the above chemical 
network, we could safely neglect it when dealing with the ionisa- 
tion radiation transfer, ionisation heating and cooling for simplicity, 
because of the much lower abun dance of the Helium compared to 
hydrogen atoms (Dvson & Will iamsll997h . The treatment of ionis- 
ing radiation from nearby stars into the above SPH code is based 
on solving the following ionisation rate and radiative transfer equa- 
tions. 



dn^ 
dt 
dj 

— = -crn(\ - x)J 

dz 



(5) 
(6) 



where n^ = xn is the electron density by ionising hydrogen atoms, 
and X is the ionisation fraction; I = o-n{\ - x)J is the ionisation 
rate, J is the flux of Lyman continum photons at the site and cr is 
the ionisation cross-section of hydrogen in the ground state. K = 
njag = x^n^ag is the recombination rate, where ag is the effective 
recombination c oefficient under the assum ption of the 'on the spot' 
approximation ( Dvson & Williams! 1 997b . In its original definition, 
the recombination coefficient 



includes all of the individual recombination coefficients a, to the 
hydrogen atomic energy level i. In the 'on the spot' assumption, 
recombinations into the ground level (/ = 1) do not lead to any 
net effect on the change in ionisation rate, since the photons re- 
leased from this recombination process are able to re-ionise other 
hydrogen atoms on the spot. Therefore ci-i can be neglected and the 
resulting net recombination coefficient can be written as 



The dependence of the re combination coefficient an on the temper- 
ature can be expressed bv lHummer & SeatonI il963h . 

as = 1.627 X IO-'^f;"-(l - 1.657/ogio4 + 0.584/^^) (7) 

where t, = 10"*r (T in K). 
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The equations |6| is numerically integraded at each time step 
and with the density distribution n[]{x,y,z) which resulted from 
the dynamical evolutions. where Jq is the flux of Lyman contin- 
uum photons at the surface the molecular cloud and we take the 
value 7o = 1.5 x 10" cm~^s~', which is the same as that used by 
IWilliams etanj200lh . 



2.2.2 lonisation heating and recombination cooling 

The implementation of the thermal model for the lonisation process 
includes the energy input by photoionisation and energy loss by 
recombination. When an atom with lonisation energy £[ is ionised 
following the absorption of one photon of frequency v and energy 
E = hv, it releases one electron which carries the excess kinetic 
energy Ey-Ei which will be transfered to the gas through collisions 
with other gas particles. The resulting heating rate is described by 
the following equation 



;j(l - x)crJkT, 



(8) 



where k is the Boltzmann constant and T, is expressed as 
jCanto et aliggj) 



xl + 4xo + 6 
x^ + 2x0 + 2 



(9) 



where T^ff is the stellar temperature, and xq = (hvo)/(kTeff) with 
Vo being the frequency of the Lyman limit. 

When a free electron in the plasma is captured by a proton, a 
photon is emitted and an amount of energy £i + m^vj /2 is removed 
from the i nternal energy of the gas . The cooling rate due to this 
process is iHummer & Seatonll963h 



A,.„. = Pb n^x^kT 



(10) 



where / 



■■ as X (1 +0.158r,). 



2.2.3 Cooling by collisionally excited line radiaiton 

Accompanying with the lonisation of hydrogen atoms, oxygen is 
ionised to (Oil) as well. Consequently, the collisional excita- 
tion of low-lying energy levels of Oil makes a significant cooling 
elfect in spite of GIFs low abundance, since Oil has energy levels 
with excitation potentials of the order of kT . We use the follow- 
ing simplified formula derived by Raea et al ( 2002) to calculate the 
cooling rate due to the collisional excitation of Oil, 



Aco/ff - A„;;;(i) -I- AcolliQ.) 

with 

tog,o[^^]= 7.9?, -26.8 

n^noii 

KoUiQ) 

log\o[ ] = 1.9: 

rierioii 



■20.5 



(11) 



(12) 



(13) 



with ?i = 1 - 2000 K/r and = 1 - 5 X 10-* K/r and no,i = xng. 
The above formula is valid under the limit of low electron density 
«j. < 10* cm"', which is true for the BRCs systems. 



2.2.4 The pressure of the cloud and the boundary conditions 

At the start of a simulation, the molecular cloud is assumed to have 
a spherical symmetry which is the general assumption adopted by 
most existing theoretical models as the initial morphology of the 




Figure 1. The initial geometric and the surrounding UV radiation config- 
urations. The heavy arrow lines represent the strong Lyman continum flux 
from the massive stars above the molecular cloud and the short arrow lines 
represent the isotropic FUV radiation (hv < 13.6 eV) flux from the inter- 
stellar medium environment. 



molecular cloud. The geometrical evolution of an initially spheri- 
cal molecular cloud under the effect of external environments may 
provide a platform for us to investigate the influence of the external 
radiation field on the morphology of the cloud. For investigation 
of the star formation at the tips of Eagle Nebula fingers, it is the 
dynamical evolution of the head structure that we are interested in, 
therefore the above assumption on the initial shape of the simulated 
molecular cloud should be reasonable. The mass of the cloud is uni- 
formly distributed through a sphere of radius R. In the simulation 
we choose the middle finger of the Eagle Nebula as an example 
since it has a simpler structure than the other fingers iWhite et alJ 
[1999). The initial geometric and UV radiation field configurations 
of the simulated cloud are shown in FigureQ 

For the head structure of the middle finger of the Eagle Neb- 
ula, with a mass ~ 31 Mq and a temperature of 20 K, the Jeans 
length Rj = 0.2pc. However the radius of the head structure is 
about 0.1 pc, which menas that the structure is unstable against the 
self-gravity of the cloud according to Jeans criteria. A simple SPH 
simulation reveals that it would collapse in < lO^years. 

However it is well known that molecular clouds are supported 
by turbulent and magnetic pressures in addition to the thermal pres- 
sure P,i, = nkT. In the case of the Eagle Nebula, the observational 
results conclude that a large-scale ordered magnetic field is not 
likely to to produce sufficient internal pressure balancing the exter- 
nal pressure and there may be a disordered component of magnetic 
field which provides an isotropic pressure to balance the external 
pressure in the head structure ( White et al. 1999.1 The relevant tur- 
bulent pressure P,i, and the induced isotropic pressure Pg due to 
magnetic field are related to the thermal pressure by two dimen- 
sionless parameters I Gorti & HollenbaclL20021 



A*. 

P,h' 



(14) 



The internal pressure of the molecular cloud then consists of 
three terms: 



P = {\+a+p)P„ 



(15) 



Observations of molecular clumps in star-forming regions 
Jjiiina etalj|l999h indicate that the values of a range from about 
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in cold, dark clouds to > 2 in regions of massive star forma- 
tion. There is also observational evidence to suggest that turbulent 
support in clumps decreases on the smaller scales of star-forming 
clumps (r^o ^ 0.2pc), where a < 2 (Goodman et al 1998). There- 
fore we adopt the value of a = 2 in our simulations considering 
the actual diameter of the condensed head at the tip of the middle 
finger is 0.2 pc. 

Measurements of magnetic fields and hence p are difficult to 
make, but present observational data sugge st a wi de range of values 
for yS, ranging from to a few ( Crutche3 ll999l) . We take j8 as a 
parameter in our simulations and determine a proper value of p for 
the Eagle Nebula so that the total internal pressure can balance the 
self-gravitation when the UV radiaition has not been switched on in 
order to observe the effect of UV radiation triggered star formation 
in a molecular cloud. A value of y8 = 6 is then obtained in this way. 

We start with a warm spherical cloud of mass 68 Mq and tem- 
perature of 60 K. The initial radius of the cloud is 0.4 pc and is 
assumed to have an uniform mass distribution and zero velocity 
distribution. The boundary condition adopted is of constant exter- 
nal pressure which corresponds to an external medium composed 
primarily of atomic hydrogen with n(HI)=10 cm"^^ and T = 100 K. 

With these implementations to the original SPH code, we have 
developed a comprehensive three dimensional SPH code of the ra- 
diation driven implosion effects of young stars on their natal molec- 
ular clouds, which allows a detailed investigation on both dynam- 
ical and thermal evolution of a molecular cloud under the effects 
of nearby massive stars in a more accurate and realistic way than 
other existing models, because our model self-consistently includes 
most of the necessary physical processes which are universal in as- 
trophysical environment. 

In the following subsections, we present and discuss the re- 
sults from the simulations based on the above described model with 
20,000 simulated particles applied. 



3 RESULTS AND DISCUSSION 

3.1 The dynamical evolution of the simulated molecular 
cloud 

Figures 1^01 and Bl together show the number density evolution of 
the spherical molecular cloud with the initial conditions stated in 
Section 3.2.3. The overall evolutionary process can be categorised 
by three different stages according to the configuration of the simu- 
lated molecular core. In order to distingush different morphologies 
of the core, we define a flat/elongated core as an object whose hor- 
izontal dimension is bigger/smaller than the vertical dimension. 

Flat core formation stage: The left panel in Figure |5| shows 
the number density distribution of the simulated cloud at an early 
time t = 4300 years. When the intensive Lyman continum flux from 
nearby massive stars falls onto the surface of the upper hemisphere 
of the cloud (front surface) as shown in Figure the gas within 
a thin layer at the surface of the upper hemisphere is ionised and 
forms an ionisation front at the front surface. The ionised gas is 
heated during the ionisation and the temperature increases. The in- 
creased pressure due to the temperature increase at this top layer 
drives an isothermal shock into the cloud, which compresses the 
neutral gas ahead of it (i.e., below the front surface). The contour 
lines at the front surface layer of the cloud in the left panel in Figure 
|5|show that a density gradient has built up there. The increased den- 
sity in the neutral gas leads to a rapid mass accumulation to form 
a condensed core. In Figure|5| the accumulated mass in a volume 



of 0.2 X 0.2 X 0.35 pc^ ( similar dimensions to the condensed core 
at the head of the middle finger of Eagle Nebula) centered at the 
densest point in the cloud, is shown as a function of time over the 
whole simulation period. The mass accumulation rate in the speci- 
fied volume is about 93 Mq My"' in the first few thousand years. 
On the other hand, the ionised and heated gas at the cloud surface 
flows radially away from the surface of the cloud and forms an 
evaporating layer surrounding the surface of the cloud, shown as 
the 'atmosphere' outside of these countour lines in the left panel of 
Figure |2| 

The above stated shock compression and gas evaporation can 
be observed from the surface of the rear hemisphere as well, but 
the effect is much weaker than that at the front hemisphere, because 
the heating effect of photoelectric emission of electrons from dust 
grains is much weaker than that of hydrogen ionisation by Lyman 
continum radiation. 

At time t = 0.666 x 10^ years, the compressed layer (shock 
front) leading the ionisation front moves toward the rear hemi- 
sphere, while a flat core starts to form in the front hemisphere of 
the cloud due to the shock, as shown in the middle panel of Fig- 
ure|2| During the next 5x10"* years, the core is continously com- 
pressed and becomes flattened due to shock compression when 
t = 1.107 X 10' years as shown in in the right panel of Figure 

HI 

With further compression of the core, as shown in the left 
panel of Figure|3| the small 'ear-like' structures (indicated by the 
white arrows) similar to those in the iLefloch & LazareS (19941) 
simulations appear from the two sides of the cloud structure where 
the compressed layer meets the gas particles at the rear edge of 
the cloud at ~ 90°. A small 'nose' structure also grows from 
the rear edge of the cloud along the central line x = 0, which is 
due to the maximum ionisation effect along the surface normal di- 
rection. Although the appearance of the 'nose' structure from the 
central line, the morphology of the cloud and the core are still in a 
flat shape which is outlined by the contour lines in the left panel of 
Figur^ In the next few 10* years, the 'ears' gradually merge with 
the 'nose' structure in the center and when t = 2.564 x 10' years, 
the 'ears' finally disappear and a flat core with a 'nose' structure 
along the center line, is shown in the middle panel in Figure^ Now 
the whole cloud is in a quasi-spherical morphology. In the section 
13.2.11 we will present a detailed discussion about the mechanism 
for the formation of the 'ear' and 'nose' structures. 

Elongated core formation stage: In the next stage of the evo- 
lution of the simulated cloud, the whole molecular cloud structure 
gradually becomes elongated, as shown in the right panel of Figure 
0when t = 3.066 x 10' years. When the time t = 3.418 x 10' years, 
as shown in the left panel of Figure a finger structure appears. A 
dense and elongated core having a radius of 0. 1 pc in the plane per- 
pendicular to the ionising flux direction and length of 0.35 pc along 
the ionising flux direction forms as shown in the left panel of Fig- 
ure|4| The mass included in the above core is about 31 Mq, and the 
average hydrogen number density is about 10' cm"^. The configu- 
ration of the simulated molecular cloud at this stage is very similar 
to the head structure of the middle finger observed bv lWhite et alj 
( 1999). We will further discuss other physical features of the cloud 
structure at this moment of time in the following sections. 

Triggered star formation stage: When t = 3.675 x 10' years, 
a highly condensed spherical core starts to form at the head of the 
elongated structure due to the self-gravitation of the core, as shown 
in the middle panel of Figure |3] With further condensation of the 
gas into the center of the spherical core, the hydrogen density n 
in the center of the core increases dramatically, reaching 10** cm"^ 
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Figure 2. The evolution of the Hydrogen number density n (cm" ) of the simulated molecular cloud from early stages t = 4300 year to ( = 0.11 Myr, when a 
flat core has formed in the front part of the cloud. The axis z is along the strong radiation flux from the massive stars, and the x axis is in the plane perpendicular 
to z. The initial center of the cloud is at (x, z)=(0.0) 
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Figure 3. The evolution of the number density n of the simulated molecular cloud over the period of ? = 0.2 - 0.3 Myr, when the core gradually becomes 
elongated. The white arrows point to the place where the 'ear-like' structures are. 



when time t = 3.753 x 10' years. The mass accumulated in the con- 
densed elongated core (0.2 x 0.2 x 0.35 pc^) reaches 34 Mq when 
t = 3.753 X lO' year through compression by the isothermal shock 
induced by the ionising radiation, which causes the mass to con- 
dense and then collapse into the center by gravitation, i.e., triggers 



the next generation of star formation in the center. The material in 
the condensed core quickly collapses into a smaller sphere of radius 
of 0.06 pc and a mass of about 25 M©. 

The final mass of the whole finger structure is 60 M©. In the 
whole evolutionary process, the molecular cloud has lost a mass of 
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Figure 4. The final evolutionary stages of the condensed core collapsing into triggered star formation after 0.375 Myr 
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Figure 5. The Evolution of mass accumulated into the core of dimension 
0.2 X 0.2 X 0.35 pc^ centered in the densest point in the simulated molecular 
cloud. 



8 Mq throungh photoevaporation over 0.39 My, which results in 
an average mass loss rate of ~ 20.5 Mq My"'. This value is 
similar to 21 Mq My~' wh ich was obtained by Lefloch's simulation 
In the next section, we will discuss the 



Vly ' wh i( 



iLefloch & Lazarefl 
mechanism for the appearances of different morphologies of the 
cloud during the evolution of the cloud. 



3.2 The dynamical evolution of the cloud structure 

By observing the changes in the morphology of the cloud structure, 
we are interested in the features of growth of the 'ears' and 'nose' 
structures and also the time durations for the cloud structures to 
stay in the 'flattened core' and 'enlongated core' stages. We will 
discuss them by examining the evolution of the configuration of 
the velocity field shown in Figure|6| which reveals an approximate 
evolutionary picture of the momentum of the cloud structure. 



3.2.1 The growth of the 'ears' and 'nose' structures 

At early evolutionary stages, while the ionised gas at the top layer 
of the front surface of the spherical cloud is heated and evaporated 
from the surface, an isothermal shock travels in the opposite direc- 
tion towards the rear part of the cloud and compresses the neutral 
gas ahead, shown as in the top left panel in Figure|5| The high pres- 
sure induced velocity distribution V-{G) at the front surface layer 
is appr oximately given by (Lefloch & LazarefF. |Lefloch & Lazarefll 



VziO) = VM(cos(e)yi' 



(16) 



where 6 is the angle at a surface point r = R from the z direction, 
and V^(0) is the velocity at the the point r = R,9 = 0, which is 7.5 
pc My"' (7.4 km s"') at the time t = 4300 years. The gas close to 
the symmetry axis (x = 0) has a smaller 6 at the front surface layer 
will move foward faster than gas which lies further away (larger 
9 ) from the axis x = 0. This velocity gradient over 9 causes the 
deviation of the morphology of the front surface of the cloud from 
the hemisphere, i.e, the front surface of the hemisphere becomes 
squashed, which can been seen from the middle panel in Figure|2| 
The gas particles close to the outermost parts of the cloud (9 ~ 
90°), as shown in the top-middle and top-right panels in Figure |S| 
obtain an extra velocity along the z-direction due to the non-radial 
isothermal shock compression, so that they intersect the rear edge 
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Figure 6. The velocity field evolution of the simulated cloud over a period time of 0.37 Myr. The length of the arrows indicates the magnitude of velocity and 
the longest arrow represents a velocitiy of +28 km s" ' . 



of the cloud and form the 'ear-like' structures on the two sides of 
the cloud, which is shown in the left panel of Figure^ 

Also the gas near the symmetry axis (x = 0) progresses for- 
ward to the rear hemisphere much faster than that further from the 
symmetry axis due to its higher velocity obtained from the radi- 
ation induced compression, so that a 'nose' gradually grows out 
from the shock front along the x = Q axis over the next 0.44x10^ 
years, dragging the most of the mass of the cloud with it, as shown 
clearly in the right-top panel of Figure|6| 

In the following evolutionary stages, after the gas in the ears 
passes the front hemisphere they will merge into the moving nose 
structure along j = axis because the effect from the ionisation 
flux is greatly decreased inside the core region and the gravitational 
force of the nose structure becomes more important so that the gas 
in the 'ears' converges into the 'nose' structure. This effect can be 
seen from the bottom-left and bottom-middle panels of Figure|6|by 
the directions of the velocities of the gas in the 'ears' pointing to the 
nose. The 'ears' merging with the moving 'nose' gradually creates 
the tail structure of the cloud complex. 



The average velocity gradient in the dense core in the left 
panel of the Figure |31 (which corresponds to the current state of 
the observed head structure of the Eagle Nebula) is |j: ~ 1.88. km 
s"' p c~' , which is consi stent with the observed value of 1.7 km s"' 
pc"' iWhite et alJl999^ . The compressed gas left in the head of the 
cloud complex now has reached a density which exceeds its Jeans 
density so that it collapses into a dense core to form a new star 
there, while the whole cloud complex keeps moving forward along 
the -z direction due to the 'Rocket effect' iOort & SpitzeJlT955h 
arising from the evaporating gas leaving the front surface of the 
cloud, as shown in the bottom-right panel of Figure|5| This 'Rocket 
effect' makes the mass center of the cloud complex move a distance 
of 0.4 pc over a period of 0.37 My. 



3.2.2 The duration of the flat and elongated morphologies 

Some of the above morp hological features can be found in the 
iLefloch & LazareS ^1994^ model as well, e.g., the 'ear' structures, 
the elongated cloud complex (or cometary ) morphology formed at 
the final stage of the evolution of the cloud. However we can not 
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ignore some profound differences in the simulated morphologies 
presented by these two models. Firstly in our simulations, the 'ear' 
structures does not last as long as in Lefloch & Lazareff's model, 
which we think may be the consequences of a) inclusion of grav- 
itation in our model and b) the 3 dimensionlity of our model, so 
the gas in the 'ears' in our model converges into the symmetry axis 
X = much earlier than that in Lefloch & Lazareff's model (which 
is clearly seen by the directions of the velocities of the gas at two 
sides), which results in the fact that the cloud complex stays longer 
at a flat-core morphology in our simulations . As can be seen from 
Figure |2j0 the ratio of the time for the cloud complex staying in 
flat-core morphology to that in the elongated-core morphology is 
~ 11 ~62%, which means that the simulated bright rimmed cloud 
spends 63% of its lifetime in a type A (with a flat core) morphol- 
ogy, and ~ ~37% of that in type B and C (with an elongated 
core), before it devolops further to form cometary structure or col- 
lapse to form a new star. We reach a similar conclusion from the 
results of a series of simulations performed for clouds of different 
masses. According to the statistical average principle, the above 
results means that we should expected to find twice as many of 
flat-core BRCs as enlongated-core BRCs when we observe a big 
group of BRCs in a similar astrop hyjical environmen t. In c ompar- 
ison, the relevant 2-D modelling <Lefloch & Lazare fi ll994ill995l) 
gave a value of ~ ~34% for the time duration for the type A 
morphology. 

The above simulated morphological sequence of the cloud 
compl ex h as been validated b y many observations. Suaitani et al. 
il99ll) and ISugitani & Qgural ( 1994) classified their observed 89 
BRCs (the SFO southern catalogue) into types of A, B and C with 
an increasing curvature of clouds' rim, shown as in the Figure 
which also roughly corresponds to the three stages in our simu- 
lated BRC evolutionary morphologies.lt is found that the majority 
of BRCs (63%) in the SFO southern catalogue are type A clouds, 
the remaining 37% in type B and C. Our latest radio continuum and 
molecular line observational study of bright-rimmed clouds reveals 
a v ery similar B RC morphological distribution to that of Sugitani et 
al iUrauharlTet al 2003, 2005.). A furthe r investigation on the struc- 
ture of BRCs. logura & Su girani|l 998|) supported bv p hotoionisa- 
tion induced shock models j Vanhala & Cameroni) 19981) . also sug- 
gested that the above three morphological BRC types are possibly 
a time evolutionary sequence. Therefore the above interpretation of 
our simulated results on the morphological evolution of a BRC can 
be well supported by observations. The value of 63% for the ratio 
of the number of flat-core BRCs to that of elongated-core BRCs is 
consistent with the observational value 66%. 

In order to further assure the role of the gravity played in 
the triggered star formation process within the BRCs, we repeat 
the simulation without including gravity in the code , the r esult 
reveals similarities to those from iLefloch & LazareS il994) and 
[williams et al. ( 2001 ) simulation, i.e., a) the core has a shorter time 
duration in a flat core morphology than that with gravitation in- 
cluded; b) the core will not collapse and after 4-5 10' years it stays 
in an equilibrium state with an average number denisty of 10' cm"^ 
and a central density up to 10''cm"^. Therefore self-gravitation of 
the cloud indeed plays an important role in triggered star formation 
under the influence of UV radiation field. We will present more 
simulation results and analysis in our next paper for the discussion 
of the relationship between the morphology of the BRCs to cloud's 
initial structures. 



3.3 The propagation of shock and the inonization front 

An important issue in the investigation of the star formation within 
the Eagle Nebula fingers is to examine whether the propogation 
of the shock induced by the ionising radiation into the structure 
has triggered an irreversible collapse of the molecular cloud into a 
condensed center. From the simulated results, it is possible to study 
the propagating process of the shock by observing the evolution of 
the number density distribution n in the cloud structure. 

In Figure|8| snapshots of the distribution of the number den- 
sity n and the ionisation ratio x along the center line of the structure 
at different times are shown. The left panel of the first row shows 
that at a very early stage (/j = 0.0043 Myr), the gas in the cloud is 
pushed towards the center not only from the front surface but also 
from the rear surface due to the surrounding weaker UV radiation 
from the interstellar medium. The ionisation ratio x at the front sur- 
face reaches a peak value of 0.58. Since the ionisation ratio x is lin- 
early related to the ionising flux Jo, a sharp ionisation front is built 
up instantly. The hydrogen ionisation heating greatly increases the 
temperature at the front surface which results in a high pressure so 
that a strong compressive shock is formed at the front surface of the 
structure, as shown in the right panel of Figure|8| The bottom two 
panels show that the shock continues leading the ionisation front 
into the cloud center at f 3 =0.11 and 14 = 0.2 Myr respectively. 

Overall, the head of the finger becomes denser and denser un- 
der the effect of the ionisation induced shock compression. With 
the increased density towards the center of the head, the ionising 
flux decreases very quickly (J(z) ~ 70^ °^"*' *'^), which results in a 
rapid drop of ionisation ratio x inside the head structure. 

If we define the region with highest density in the density 
profile as the shock front, the average speeds of the shock front 
are v, = 0.33, 0.22 and 0.15 km s"' over the three periods of 
h - hJi - h and - ?3, while the corresponding speeds for the 
ionisation front are v, = 0.18, 0.165 and 0.056 km s"' respectively. 
It is shown that the propagations of both shock and ionisation fronts 
slow down because of the increased density ahead. Since the ioni- 
sation front speed i', < v,, so the gap between the ionisation front 
and the shock front slowly increases. After / ~ 0.36 Myr, the gap 
reaches a value of 0.125 pc, which is consistent with the fact that 
the peak submm continum emission from the clumps is located at 
about 0.1 pc deeper into the fingers than their photoionised front 
surface IW hite et alj ri999). In the following section, more details 
about this displacement will be discussed. 



3.4 The evaporation of the ionised gases - HII region 

While the ionising radiation induced shock precedes the ionising 
front into the cloud structure, the ionised hydrogen atoms at the 
top layer of the upper-hemisphere are heated, evaporated and move 
away from the front surface of the cloud to form a HII region 
around the front surface of the cloud. Although the gas at the other 
sites of the surface other than the front sites are heated as well by 
the surrounding environmental FUV radiation through the photo- 
electric emision from the surface of dust grains, the heating effect 
is much weaker (due to the low abundance of the dust grains) when 
compared with that of hydrogen ionisation so that the evaporating 
envelope is much thinner than that surrounding the front surface, as 
shown in Figure l2P1 

The panels in the upper row in Figure |9| show the movement 
of this evaporating layer for the first few snapshots of / = 0.0043, 
0.009, and 0.015 Myr. The average evaporation speed is about 25 
km s"' . From the bottom panels of the same figure, we can see that 
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Figure 7. The schematic of three BRC morphologies. The vertical arrows represent the incoming ionising photon flux and horizontal arrow expresses the 
direction of the possible time-evolution of the cloud complex under the influence of the ionising radiation flux. The three types A, B and C correspond to the 
three morphologies in three different stages over the evolutionary of the cloud complex. 
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Figure 8. A shock (soUd Une) is leading the ionisation front (dot-dashed Une) moving into the molecular cloud. From left to right, up to the bottom are the 
density distributions in the cloud along the synmietrical axis x = for the time t=0.0043, 0.067, 0.11, 0.2 Myr respectively. 
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Figure 9. Top row: The ionised hot gas is leaving the upper semi-spherical surface. The contours are the values of the ionisation rate x, the biggest value is 
0.58. The left, middle and right panels are the ionisation rate distributions in the evaporating layer over the head of middle finger of the Eagle Nebula at times 
0.0043, 0.009 and 0.015 Myr respectively. Bottom row: The temperature (Log T (K)) distribution in the ionised hot gas over the same period of times. 
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Figure 10. The temperature distributions at the initail and late stages of 
the evolution. Dashed line is for t = 0.0044 Myr, and the solid line is for 
/ = 0.32 Myr 

the temperat ure of this evaporati ng layer is about several 6 x 10^ 
K, fo r which Heste r et alj (l99g) estimated a value of lO"* K, and 
IWhite et aL (il999i) gave a value of 6700 K. The characteristics of 
the modelled HIT region surrounding the front surface of the Eagle 
Nubula finger fits well with the observed values. 



3.5 Evolution of the temperature profile 

The evolution of the temperature distribution of the head structure 
at times / = 0.0044 (the dashed line) and 0.34 Myr (the solid line) 
along the center line X = 0, are shown in Figure fTol It is seen 
that shortly after the ionising radiation is switched on, i.e., when 
t = 0.0043 Myr, the top layer at the front surface is heated as a 
consequence of hydrogen ionisation heating, so that the tempera- 
ture increases to ~ 6000 K. Inside the structure the temperature is 
about 60 K, close to the inital temperature of the cloud, because the 
elfect of the radiation has not yet penetrated. After the cloud has 
evolved for / = 0.32 Myr, the temperature in the surface layer of 
the cloud is about 200 K, since the ionisation flux has greatly de- 
creased at the surface of the cloud due to the ionisation absorptions 
in the surrouding HII regions (which is moving away from the main 
structure) for the recombinations of the electrons with the hydro- 
gen ions. The temperature at the rear surface increases to a similar 
value, due to the FUV radiation from the surrounding interstellar 
medium. The combination of the value of the surface temperature ( 
~ 200 K) and the hydrogen density (n ~ 5 x 10"*) yields a thermal 
pressure in the surface layer of the = ~ 10' cm"^ K. The 
total pressure j ~ 7 x 10' cm"' K according t o Ea.<15|. wh ich is 
in a good agreement with the observed value bv lPounjil998l) . The 
temperature in the core region of xxy x z = 0.2 x 0.2 x 0.35 pc', 
which is centered at (x,y,z) = (0,0, -0.43) pc, decreases to 13-18 
K, below its original value of T = 60 K, due to the enhanced radi- 
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ation extinction, and the increased cooling efficiency in this higher 
density region. 

In summary, the temperature distribution at the time t = 0.34 
Myr describes a low temperature (13 -18 K) core surrounded by 
an outer shell of warm material of temperature of 200-250 K, 
which matches well with the the temperature characteristics of the 
oberved head structure at the m iddle fingertip in the Eagle Nebula 
jWhite et alJi99^|PounJl99l . 



3.6 The main chemical abundance evolutions 

The initial conditions for all of the chemicals abundances in- 
cluded are the same as those described in the previous paper, with 
X(C+)=A'(CO)=X(HCO+)=0 and X (C)= X(0) = 10' ^ when the 
simulation starts, where Xj = -j^^ jWhite et ai]ll999h . In the left 
panel of Figure ^3 the distributions of the fractional aboundances 
of carbon bearing species are displayed for the time t = 0.32 Myr 
In the densest region between z=- 0.5 and - 0.3 pc, carbon is mainly 
in the form of CO, which correlates with the density of H2 because 
its production rate scale is oc «^(H2), which makes the peak deep 
inside the structure. Outside the core z > -0.3 pc, CI and C+ are of 
comparable abundance at the front edge of the cloud strucure, due 
to an increased ionisation of CI into C^ which results in a dramatic 
increasing aboundance of the latter. The peak of the ionisation ratio 
X at z ~ -2.7 pc reflects the position of the bright rim. This is to 
say that the ionisation peak is ahead of the CO peak by a distance 
of 1.25 pc, which provides a very consistent description of the 
relationship of the CO and Ha distributions inside the cloud head 
with the obs erved displacemen t of 0.1 Ipc between CO line and H„ 
line profiles ^White et alll999h . 

The distribution of the HCO* abundance shows a non- 
monotonic variation over the region because its production rate is 
dependent on the ionisation flux, the H bearing species and CO 
species. Near the front edge of the structure at z = -0.3 pc, these 
three factors makes a peak of the HCO^ production rate. Deeper 
inside the structure, the high abundance of CO is dominant in the 
production rate of HCO^, therefore similar pattern to that of CO, 
which is consistent with the observations o f|wMteet^jl99i). 

White et al. also found that the ratio of CI/CO abundance in- 
creases in the low column density regions. The right panel of the 
figure fm shows the corresponding ratio increases dramatically in 
the front region of the head structure, where the density of the cloud 
is very low. 



4 COMPARISON WITH THE RESULTS FROM THE 
PREVIOUS MODEL 

Firstly, we would like to point out that the estimated cloud 
thermal and chemical structures presented in the previous paper 
was based on a static model under the assump- 
tion that the cloud was under hydrostatic equilibrium, i.e., the den- 
sity distribution was fixed in time. Then the temperature and chem- 
icals profiles were estimated according to the thermal and chemical 
models described in the previous SPH code, which of course were 
not actually describe the dynamical evolution of the cloud structure. 

Secondly, a dynamical simulation based the previous SPH 
code in the absence of the magnetic and turbulent pressures sup- 
port, indicates a quick collapse of an spherical molecular cloud of 



initial mass 68 Mq and radius of 0.4 pc (or that of 3 1 M(T) and ra- 
dius o f O.lpc, as we have used in our previous paper iWhite et alJ 
'199^ in much less than 10 years which is not a reasonable result 
for the Eagle Nebula complex. 

Thirdly, as we have mentioned in the previous sections, the 
previous SPH code didn't include the Lyman continum in the ra- 
diation field, and the consequent hydrogen inoisation heating and 
cooling processes so it greatly decreased the effect of radiative im- 
plosion on the bright rimmed cloud, because the temperature near 
the bright rim of the cloud can only be increased to 190 K by the 
photoelectric effect of the FUV radiation. With the inclusion of Ly- 
man continum radiation, the ionising hearting can heat the front 
layer of the cloud up to 7000 K in the first few hundred years, a 
reasonable temperature for HII region. 

Therefore the simulation based on the previous SPH code can 
only simulate the effect of photoelectric heating of FUV radiation 
on the evolutionary of the molecular cloud, not the whole evolution 
picture of how a ubiqutous molecular cloud under the effect nearby 
massive stars can evolve into a finger-like BRC with its life ending 
at either triggered star formation, formation of stationary cometary 
structure, or being evaporated totaly by the nearby massive stars' 
radiative field. As we have shown in this paper, these effects can 
only be modelled via the comprehensively extended model. 



5 CONCLUSION 

With the inclusion of the Lyman continum in the external radiation 
flux and hydrogen ionisation heating effect into the energy equation 
in the previous code, and turbulent and magnetic pressures into the 
internal pressure of molecular cloud, the upgraded 3D SPH model 
successfully describes the dynamical evolution of an initially spher- 
ical uniform molecular cloud into cometary/triggered star forma- 
tion under the influence of the strong UV radiation from nearby 
stars plus an isotropic surrounding interstellar medium radiation, 
and self-gravitation of the molecular cloud. 

The application of the model to an initially spherical molecu- 
lar cloud of mass 68 MO and a temperature of 60 K results in an 
Eagle Nebula finger-like structure with a very dense core of a di- 
mension of 0.2 X 0.2 X 0.35 pc^^ and a mass of 31 MO at its head 
after 0.34 Myr. The core has an average density of 10' cm"^^ and 
a temperature of about 13-18 K with a velocity gradient of 1.91 
km s"' at time t = 0.34 Myr. The dense core is surrounded by 
a warmer layer of temperature 200-250 K. The simulated HII re- 
gion is of temperature betweem 3-7000 K. The charateri sties of 
the simulated chemicals distributions are: a) CO abundance peak is 
0.125 pc deeper into the core from the peak of ionisation fractional 
ratio X at the edge of the head structure; b)HCO+ abundance fol- 
lows a simillar pattern to CO in the core region; c)the abundance 
ratio of CI/CO increases with the decreaes of the cloumn density 
of the cloud. The above described characteristics of the simulated 
head structure coincide very well with the recent observational data 
(White et al. 1999). The simulation results also shows that the core 
will collapse into its center after another 7800 years. 

It is concluded that the radiation induced compression wave 
has yet to pass the head structure and that the currently observed 
head structure at the tip of the middle finger is a part of an ob- 
ject which evolved from an initially larger and warmer pre-existing 
molecular clump. This cloud clump has been under the effect of a 
strong UV radiation from nearby stars for about 0.34 Myr and is 
currently at a transition stage towards triggered star formation. 

The simulation reveals some of the similar morphological 
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Figure 11. The main chemicals distributions at the time / = 0.34 Myr 



structures in the evolutionary process of the molecular cloud to 
those displayed by Lefloch & Lazareff's 2-D model and obtains 
a very similar estimation on the photoevaporation rate (20 MO Myr 
^' vs 21 Mo Myr ) from the surface of the front side, which 
shows the reliability of our modelling. This also shows that the 
Lefloch & Lazereff's 2-D model is appropriate to consider the pho- 
toevaporation of a bright-rimmed cloud, but not the issue of trig- 
gered star formation in a the cloud. 

The estimation on the relative lifetimes of between type A 
and type B/C BRCs shows that a BRC remains in type A mor- 
phology nearly twice the time spent in the type B/C morphology, 
which means that we should observe double the number of BRCs 
in Type A morpho logy. This pred i ction agree s well with the obser- 
vational results bv^ Sugitani et alj (1991") and'Urguh art et all i2003l 
(2005), thus provides a reasonable explanation for the observed un- 
balanced numbers of flat-core and elongated-core BRCs. 

Comparing the above results with those based the previous 
SPH model, we can conclude that the hydrogen ionisation heating 
by Lyman continum plays a much more important role in the dy- 
namical evolution of a molecular cloud than that of photoelectric 
heating by the surrounding interstellar medium radiation. Also the 
inclusion of the turbulent and magetic field pressures in the molec- 
ular cloud in the upgraded model makes the model more realis- 
tic so that it can trace the origin of an observed cloud, reveal the 
whole dynamical evolution process of a molecular cloud to trig- 
gered star/cometary formation under the influence of nearby stars, 
and predict the future evolution of the observed object. 
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